A Concept for A Dark Matter Detector Using Liquid Helium-4 
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Direct searches for light dark matter particles (mass < 10 GeV) are especially challenging because 
of the low energies transferred in elastic scattering to typical heavy nuclear targets. We investigate 
the possibility of using liquid Helium-4 as a target material, taking advantage of the favorable 
kinematic matching of the Helium nucleus to light dark matter particles. Monte Carlo simulations 
are performed to calculate the charge, scintillation, and triplet helium molecule signals produced 
by recoil He ions, for a variety of energies and electric fields. We show that excellent background 
rejection can be achieved based on the ratios between different signal channels. We also present 
some concepts for a liquid-helium-based dark matter detector. Key to the proposed approach is the 
use of a large electric field to extract electrons from the event site, and the amplification of this 
charge signal, through proportional scintillation, liquid electroluminescence, or roton emission. The 
sensitivity of the proposed detector to light dark matter particles is estimated for various electric 
fields and light collection efficiencies. 

PACS numbers: 34.50. Gb, 33.50.-j, 82.20.Pm 
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INTRODUCTION 



Dark matter, while evident on multiple astronomical 
length scales through its gravitational effects, has an un- 
known intrinsic nature. Data from primordial nucleosyn- 
thesis [l[ , the cosmic microwave background , structure 
formation Q, and microlcnsing observations Q imply 
that the dark matter cannot be composed of baryons 
or active neutrinos, implying new physics beyond the 
Standard Model. Experimental direct detection of dark 
matter particles, illuminating their mass and interaction 
properties, would therefore create crucial new scientific 
understanding in both astrophysics and particle physics. 

A particularly compelling model for dark matter is that 
it consists of Weakly Interacting Massive Particles, or 
WIMPs [1, Q , with the feature that a massive particle 
in the early universe interacting through a weak-scale 
cross-section yields a thermal relic abundance approxi- 
mately that observed for dark matter. Over the past 
few decades, models of WIMP dark matter have cen- 
tered on constrained minimal supcrsymmetry (CMSSM) 
models which predict a stable neutralino with mass 
greater than 40 GeV, limited to higher masses by the 
requisite mass difference between the chargino and neu- 
tralino. Also, it is commonly argued that in the con- 
text of supersymmetry it is most natural for the dark 
matter mass to be comparable to the weak scale Q . 
As a result, most direct dark matter experiments have 
been designed to have excellent sensitivity to dark mat- 
ter particles with mass comparable to or greater than the 



weak scal e, y et most of these, including the CDMS [l^ , 
ZEPLIN and XENON programs, see no ev- 

idence for such high mass dark matter particles, down 
to the recent XENONIOO spin-independent cross-section 
limit of about 2 x lO'^^ cm^ at 55 GeV [ll|. At the same 
time, the DAMA [H, CoGeNT [H, and CRESST ^ 
experiments have seen event rate anomalies that can be 
interpreted in terms of direct detection of light WIMPs, 
and a number of astrophysical anomalies may be inter- 
preted in terms of light WIMP annihilation [1^. Mean- 
while, many new theories of light WIMPs have been de- 
veloped, and this is currently an area of active develop- 
ment in particle phenomenology. Models for light dark 
matter often involve a new mediator particle as well as 
the dark matter itself, and include the next to minimal 
supersymmetic model (NMSSM) fl9|. a symmetric dark 
matter [2^, WIMP less dark matter l2 l[ , singlet scalars 
[2^, dark sectors with kinetic mixing j23l|. mirror mat- 
ter [2J]. These models can all evade constraints on light 
WIMPs from the cosmic microwave background [25l , the 
Large Hadron Collider [U 
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and Fermi-LAT 

Considerable excitement has been generated over the 
possibility that dark matter particles are relatively low 
in mass. The difficulty is detecting them, since lighter 
WIMPs have less kinetic energy and only deposit a small 
fraction of it when elastically scattering with standard 
heavy targets like germanium and xenon. 

In general it is difficult for heavy targets to be sensitive 
to light WIMPs, since for typical energy thresholds they 
are only sensitive to a small part of the WIMP veloc- 
ity distribution. Models of the WIMP velocity distribu- 
tion typically assume a Maxwcllian distribution of f{v) = 

exp — (w + veY /vq , where ve — 244 km/s is the velocity 
of the Earth around the Milky Way, and vq ~ 230 km/s 
is the virialized velocity of the average particle that is 
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gravitationally bound to the Milky Way [28| . This distri- 
bution is expected to be roughly valid up to the Galactic 
escape velocity v^sc — 544 km/s, above which the veloc- 
ity distribution is zero. A plausible energy threshold for 
Xe, Ge, and He dark matter experiments is about 5 keVr. 
But for a 5 GeV WIMP, such as predicted by asymmetric 
dark matter models j20l | , its velocity must be particularly 
large to deposit at least 5 kcV. This minimum velocity, 

Vmin, is equal to Vrain = \J\- ^R- Mj- j T , whcrc Er is 
the recoil energy, r is the WIMP-target reduced mass 
r = Md-Mt/{Md + Mt), Md and Mr are the masses of 
the dark matter particle and the mass of the target nu- 
cleus, respectively. For Er of 5 keV and Md = 5 GeV, 
Vmin is equal to 1127, 864, and 427 km/s for Xe, Ge, and 
He respectively. So for this example, Wmm is above Uesc 
for Xe and Ge, but not for He. The lower limit of the 
WIMP-target reduced mass that a detector is sensitive 
to is given by 
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runvit = ^/EfMT/2, 



(1) 



where Et is the energy threshold. So a kinematic figure 
of merit for light WIMP detection is the product of the 
energy threshold and the target mass, which should be 
minimized for the best light WIMP sensitivity. 

This challenge of combined low energy threshold and 
low target mass can likely be met through the use of 
liquid helium as a target material. In this paper we in- 
vestigate the use of liquid helium as a target for light 
dark matter particles in the mass range of 1 to 10 GeV. 
In Section II we outline the properties of liquid helium 
in the context of particle detection, in Section III we de- 
scribe possible configurations of helium-based detectors. 
The detector can be operated at T '--^ 3 K, adopting pro- 
portional scintillation or electroluminescence for charge 
readout; or it can be operated at T ~ 100 mK using 
bolometers for light and charge readout. In Section IV 
we examine the sensitivity of liquid helium detectors to 
hght WIMPs. We conclude in Section V. 



II. LIQUID HELIUM AS A DETECTOR 
MATERIAL 



Superfluid helium has been used for a detector mate- 
rial for many applications. Most detector concepts take 
advantage of the special excitations of the superfluid, 
and involve detection of phonons, rotons, or quantum 
turbulence. One example is the HERON concept [2^ 
for pp-solar neutrino detection with rotons in superfluid 
helium-4 at a temperature of -^100 mK. The HERON 
researchers also considered using such an instrument to 
look for dark matter [s^, iSlJ i with the possibility that the 
roton/vortex generation by electrons in an applied elec- 
tric field, combined with prompt roton detection, could 
be used for particle discrimination. Also, the roton sig- 
nal could carry information about the nuclear recoil di- 
rection. Another is the ULTIMA concept 32| for dark 



matter detection with quantized turbulence in superfluid 
hehum-3. Both of these concepts have been the subject 
of considerable research and development in the past few 
decades. 

Along with its many unusual properties related to su- 
perfluidity, liquid helium also produces substantial scin- 
tillation light and charge when exposed to ionizing radi- 
ation, just like liquid xenon and liquid argon which are 
already used extensively in the search for dark matter. 
Some ultracold neutron experiments already make use 
of the prompt scintillation of liquid helium; for exam- 
ple the measurement of the neutron beta-decay lifetime 
[33| and search for the neutron permanent electric dipole 
moment at the Spallation Neutron Source [ll, HI] . The 
prompt scintillation yield in liquid helium is well known, 
measured by the HERON collaboration to be about 20 
photons/keV electron equivalent (keVee). 

Depending on particle species, energetic particles elas- 
tic scattering in helium can lead to electronic recoils 
(gamma ray, beta scattering events) or nuclear recoils 
(neutron or WIMP dark matter scattering events). The 
recoil electrons or He nucleus collide with helium atoms, 
producing ionization and excitation of helium atoms 
along their paths. The ionized electrons can be extracted 
by an applied electric field. The decay of the helium ex- 
cimcrs gives rise to scintillation light. A fraction of the 
deposited energy is converted into low-energy elementary 
excitations of the helium, i.e., phonons and rotons. Sig- 
nals from all these different channels may in principle be 
used to detect and identify the scattering events. The 
key for dark matter detection is to be able to suppress 
the electronic recoils that make up most of the back- 
grounds from the nuclear recoils that would make up 
a WIMP signal by use of event discrimination. In this 
section, we shall estimate the nuclear and electronic re- 
coil signals due to ionization charge, prompt scintillation 
light, metastable Hcj molecules. We shall present re- 
sults of Monte Carlo simulations showing that excellent 
background rejection can be achieved for the purpose of 
WIMP dark matter detection, based on the ratio between 
these different signals. 



A. Low energy nuclear recoils in helium 

1. Charge exchange processes 

A WIMP dark matter scattering event in liquid he- 
lium would result a recoil helium nuclei. Depending on 
the energy involved in the scattering process, the recoil 
He can be a bare ion (He^"*") or a dressed ion (He^+), 
or even a neutral helium atom (He°). The recoil He dis- 
sipates its kinetic energy through collisions with ground 
state He atoms. Such collisions can be elastic or inelastic 
that lead to ionization or excitation of He atoms. The 
ionization and excitation cross-sections are different for 
the recoil He ion in different charge states. As the fast 
recoil He ion slows down, interactions involving electron 
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FIG. 1: (color online). Charge exchange cross sections due to He'', He'^"'", and He^^ interacting with ground state He atoms. 
The curves were fitted to experimental data by polynomial functions, (a) cr'^^: — (this work), • (Refjs^), o (Ref. 131), 



▲ (ReM3M|); a"^: (this work), T (Refjsdl). (b) — (this work), ▲ (Ref. [Sg), □ (Ref.j40|), o (Ref. 

(Ref. a^': (this work), T (Ref. [ij), x (Ref. [U). (c) a'^°: — (this work), ▲ (Ref. [U), □ (Ref. 

(Ref. H), * (Ref. [13); o-^^: (this work), • (Ref. H), (Ref. [!§]),♦ (Ref. H), A (Ref. [H) 
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capture and loss by the projectile become an increasingly 
important component of the energy loss process. Charge 
transfer can produce residual ions without the release 
of free electrons, and free electrons can be ejected from 
the moving ion (or neutral) with no residual ions being 
formed. 

Charge transfer cross sections are generally designated 
as a'^^ where i represents the initial charge state of the 
moving ion, and / is the charge state after the col- 
lision. For a complete description of the full slowing 
down of a recoil He, we need cross sections for one- 
electron capture cr^^ and two-electron capture (t^° for 
He^"*" , one-electron capture cr^° and one-electron loss cr^^ 
for He^+j and one-electron loss (t°^ and two-electron loss 
cr°^ for He°. In Fig. [U we show the six charge exchange 
cross-sections based on available experimental data for 
He", He^"'' and He^"*". These cross sections were least- 
squares fitted by simple polynomial functions of the form 
log(cr*/) = J2„ C„(log £;)", where the C„'s are the fit- 
ting parameters, and E is the particle energy in kcV. 
Smooth extrapolation was carried out where the experi- 
mental data were lacking. Following the method by Al- 
lison [sO], the fractions Fq, Fi, and F2 that the moving 
particle to be found in charge state 0, 1, and 2 are given 
by 

dFo/dz = N[-Fo{a°^ + a°^) + F^a^" + Fsa^O] 
dFi/dz = N[-Fi{a^° + cri2) + Foa^i + F2(7^^] (2) 
F2 = 1 — Fo — Fi 

where TV ~ 2.2 x 10^^ cm""^ is the number density of 
liquid helium and z is the path length along the particle 
track. If the charge exchange cross-sections cr*^ do not 
vary as the He ion moves, the equilibrium charge fractions 
F^, and F^ as z ^- 00 arc given by Alhsion [H^ 



as follows: 

F^ = {ba^'^ -ga^')/D (3) 
F^ = [{a - b)a'" + gia + a'') - f{b + a^^)]/D 

in which 

a=-(aio+ai2 + a2i), b = a"^^ - 

f = - ^20, g = -(cjOi + + (4) 

D = ag-bf 
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FIG. 2: (color online) Equilibrium fractions of the charge 
states of an energetic helium ion in liquid helium. 

Fig. [2] shows the calculated equilibrium charge frac- 
tions as a function of helium ion energy based on Eq. [3] 
and Eq. S) At energy higher than a few thousands of 
keV, the helium ion appears primarily as a bare ion He-+ , 
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whereas in low energy regime (< 100 keV) the fraction 
of charge zero state He" dominates. These results are 
derived based on the assumption that cr*-^ does not vary 
as the He ion moves. In reality, since the charge ex- 
change cross-sections depend on particle energy, as a He 
ion slows down in liquid helium, the a^-^ in Eq. [2] should 
change as z varies. In this situation, a full description of 
variation of the charge fractions Fq, Fi, and F2 is given 
by the following equations 



dFo(E) _ N 



dE S{E) 
dFi{E) _ N 

dE ^ SiE) 
F2{E) = 1 - Fo{E) - FiiE) 



[-Fi(<t1 



(5) 



where S{E)=dE / dz is the total stopping power of a He 
ion in liquid helium that describes the average energy 
loss of the He ion per unit path length. S{E) is the sum 
of the electronic stopping power Se{E) (energy loss due 
to the inelastic collisions between bound electrons in the 
medium and the ion) and the nuclear stopping power 
Sn{E) (energy loss due to the elastic collisions between 
the helium atoms and the ion). Fig. [3] shows the stop- 
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FIG. 3: (color online) Stopping power of a He ion in liquid 
helium. Data are drawn from the National Institute of Stan- 
dards and Technology (NIST) database [H]. 

ping power data drawn from the National Institute of 
Standards and Technology (NIST) database [s^ • Know- 
ing the stopping power S(E) and the charge exchange 
cross-sections a^^{E), one can integrate Eq. [5]to calcu- 
late the energy dependence of the fractions of different 
charge states with a given initial condition. An exam- 
ple is shown in Fig. |4l We see that if we start with a 
bare ion He^"'' {F2 = 1) at an initial kinetic energy of 
50 keV, as the ion slows down the fractions of the dif- 
ferent charge states Fq, Fi, and F2 quickly evolve to the 
equilibrium values. This is because that due to the rela- 
tively large charge exchange cross-sections and the high 
helium number density, many charge exchange collisions 
can take place in a short path-length of the fast He ion. 
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FIG. 4: (color online) Fractions of the charge states of an 
energetic helium ion as it slows down in liquid helium. The ion 
started as a 116^"*" with initial energy of 50 keV, as indicated 
by the red circle. The arrows show how the fractions evolve 
as the particle loses its energy. 



To achieve the equilibrium charge fractions, only a few 
charge exchange collisions are needed and the energy loss 
in this process is small. As a consequence, we can safely 
use the equilibrium fractions of the charge states to study 
the slowing down of a fast He ion in liquid helium, with 
no need to consider the initial charge states. 



2. Ionization and excitation yields 



The ionization and excitation yields due to a recoil 
helium nuclei moving in liquid helium are important 
premise parameters needed for the design of a helium- 
based dark matter detector. Sato et al. [53j have studied 
the ionization and excitation yields of an alpha parti- 
cle (He^"*") in liquid helium using the collision cross sec- 
tions derived with the binary encounter theory [s^]. In 
their analysis, the charge exchange collisions are ignored 
and the fraction of the alpha particle energy that is lost 
to elastic collisions with surrounding He atoms (nuclear 
stopping) is not included. Nuclear stopping can become 
dominant when the alpha particle energy is small, which 
is known as the Lindhard effect [s^ . The energy of a re- 
coil helium nuclei in a WIMP scattering event is expected 
to be relatively low (< 100 keV). To obtain more reliable 
estimation of the ionization and excitation yields from a 
recoil helium nuclei, we present an analysis that system- 
atically accounts for both the charge exchange processes 
and the Lindhard effect. 

Assuming a continuous slowing down, the total number 
of free electrons iVei produced along the path of a recoil 
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FIG. 5: (color online). Ionization cross-sections due to He", He^"*", and He^^ interacting with ground state He atoms. The 
curves were fitted to experimental data by polynomial functions, (a) crfo^: — (this work), • (Ref. [HH]), A (Ref. [57l|). o 
(Ref. Hi), (b) al+„: — (this work), • (Ref. 0), A (Ref. [H). (c) cr^+„: — (this work), o (Ref. [HI), A (Ref. M)- 



He nuclei with an initial kinetic energy E is given by 



iVe; 



^ NdE' 
® NdE' 



S{E'] 
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ion 



Here 



and 



^Exc 



(6) 



arc the number of electrons pro- 
duced in direct ionization and in charge exchange pro- 
cesses due to He ion impact, and are given by the first 
and the second integral terms on the right side of the 
equation. (J±onj '^lom '^lon ^-re the direct ionization 
cross sections due to He°, He^+, and He^+ interacting 
with ground state He atoms, respectively. N^l'^ is the 
number of ionizations produced by secondary electrons 
that have energy higher than the ionization threshold of 
a He atom (24.6 eV). (i = 0, 1, 2) is the equilibrium 
fraction of charge state i as given by Eq. [H The ratio 
of N^i'^ to A^ei decreases with decreasing E and is only 
a few percent when i? ^ 100 keV [11]. We shall neglect 
N^i'^ in the following analysis for simplicity. To estimate 
the ionization yield, defined as Fei = Nei/E, the val- 
ues of the direct ionization cross sections are needed. In 
Fig. O the experimental data for cr?^^, a}^^, and crj^ are 
shown. We again fit the experimental data by simple 
polynomial functions log((Tion) = J2n^ni^°S -^)": ^^'^ 
extrapolate the curves where the experimental data were 
lacking. From Fig. [2] one can see that at i? < 100 keV, 
the fraction of the charge zero state (Hc°) dominates. 
The available ionization and charge exchange cross sec- 
tion data for He'' in the energy range of 0.1~100 keV 
allow us to make reliable fit and extrapolation for analyz- 
ing the ionization yield. The calculated ionization yield 
1^1 of a recoil He ion as a function of the ion energy is 
shown in Fig. [5] as the black solid curve. 

The total number of excitations iVex produced by a 
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FIG. 6: (color online) Ionization and excitation yields of a 
recoil He ion in liquid helium as a function of the He ion 
energy. 



recoil He nuclei with an initial kinetic energy E is given 

by 



. S(F)l^""-<^>-- + ^-""<^>" (7) 
f2~(-E'KVl + W.. 



where cr" , <tI^, and (7?+ are the total excitation cross 
sections due to He°, He^+, and He^+ interacting with 
ground state He atoms, respectively. Here A^ex is the 
number of excitations produced by secondary electrons, 
which can again be neglected at E < 100 keV [H^]. Ex- 
perimental excitation cross section data are limited. For 
instance, Kempter et al. estimated the excitation cross 
sections due He atom impact, but only with collision 
energy below 600 eV [62|; De Heer and Van Den Bos 
measured the excitation cross sections for He^"*" incident 
on He, but only for excitations to states with principle 
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quantum number n ^ 3 [63|. Instead of fitting the data 
to obtain the excitation cross sections, we estimate the 
excitation yield Fex = N^x/E based on the known elec- 
tronic stopping power as follows. The electronic stopping 
power Se{E) can be written as 

§ = FS'Kn (Qne + so) + criQex + {t'" + 2a°'){Qu. + XE)] 
+ FrWiin (Qh, + ei) + ^i+Q,, + ct1=(Q„, + XE)] (8) 

Here Qne = 24.6 eV is the ionization energy of a he- 
lium atom, e is the average kinetic energy of secondary 
electrons by He ion impact. Qex = J2Qij'^ij 
the mean excitation energy where Qij ^-nd aij are the 
He(i— excitation energy and the associated cross sec- 
tion, respectively. Lack of detailed information, here we 
assume Qex to be the same for the incident He ion in dif- 
ferent charge states. A = rrie/mue — 1-36 x 10"'' where 
TOg and rrins arc the masses of an electron and a He atom, 
respectively. In Eq. |8l the energy transfer model is as- 
sumed such that in a charge-loss collision, a stripped elec- 
tron is ejected from the projectile with nearly the same 
velocity as the projectile. Indeed the stripped electrons 
are observed in the spectrum of secondary electrons pro- 
duced when He ion impacts on water vapor as a peak 
centered at XE [gJ]. An energy deposition of Qne+XE 
is thus made when an electron is lost from the projec- 
tile [gHI- In an electron capture process, energy depo- 
sition is essentially due to the recoil of the ionized He 
atom and is negligible. As a result, the terms in the 
square brackets in Eq. [7] can be derived based on Eq. |8] 

= - i^o^kicn (Qne + so) + {a"' + 2a°2)(QHe + XE)] 



- FrWl+ (Q„e + ei) + T'HQne + XE)] 
-Fr['jf+iQns + e2)]} 



(9) 



Plugging Eq. [S] back into Eq. [3 the excitation yield can 
be derived as 



eJ„ s{e')q,x 

in which L is the Lindhard factor, defined as 



(10) 



1 /-^ Se{E')dE' 

EJo S{E') ■ 



(U) 



Lindhard factor designates the ratio of the energy given 
to the electronic collisions to the total energy. A plot 
of the Lindhard factor as a function of the recoil He ion 
energy is shown in Fig. [71 Since only the part of energy 
given to the electronic collisions can be used as ionization 
and scintillation signals, the Lindhard factor L is impor- 
tant for the determination of the sensitivity of WIMP 
detectors. 
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FIG. 7: Calculated Lindhard factor for a recoil He ion in 
liquid helium as a function of the He ion energy. 



In order to calculate Ygx using Eq. [TOl we need to make 
further approximations on Qex and e. Since the domi- 
nant excitation process in low energy collisions between 
He atoms and the projectile is He(ls^— !>ls2p) with an 
excitation energy of 21.2 eV [gl], we take Qex — 21 eV 
for simplicity. The average energy e of the secondary 
electrons can be expanded in power series of E. To the 
lowest order in E, we may write e ~ j{E — 24.6 eV) for 
E > 24.6 eV. Linear dependence of e on is evidenced 
for secondary electrons ejected by helium ion impact on 
water vapor with energy E < 100 keV [65|. Further- 
more, at small E, e is similar for He ion impact in dif- 
ferent charge states. We choose 7 = 0.3 for all charge 
states such that the ratio between the calculated ioniza- 
tion yield Y^i and excitation yield Ygx agrees with Sato 
et aVs result a.t E ^ 100 keV where the Lindhard effect 
is mild. Note that variation of 7 does not affect Fax at 
small E. The calculated Ygx is shown in Fig. [HI as the red 
dashed curve. 

The drop of both the ionization yield and the excita- 
tion yield at energies lower than about 50 keV is due to 
the drop of the electronic collision cross sections in this 
energy regime, as well as the loss of the He ion energy 
to elastic nuclear collisions (Lindhard effect). As a com- 
parison, for an energetic electron moving in LHe, Sato et 
al [H, lill estimated that the total ionization yield and 
excitation yield are nearly constant (Y^i^ ~ 22.7 keV~^ 

and yir'' — 10.2 kcV~^) in the energy range from a few 
hundred keV down to about 1 keV. 



B. Signals in liquid helium 

1. Charge signal 

Electrons and helium ions arc produced along the track 
of an energetic particle as a consequence of ionization or 
charge- exchange collisions. Beside these processes, ex- 
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cited helium atoms produced by the projectile with prin- 
cipal quantum number n > 3 can autoionize in liquid 
helium by the Hornbeck-Molnar process [g^I 



He* 



He 



He.^ 



(12) 



since the 2 cV binding energy of Hcj is greater than the 
energy to ionize a He(?? > 3) atom. Based on the os- 
cillator strengths for the transitions between the ground 
state and the various excited states of helium , slightly 
more than one third of the atoms promoted to excited 
states will have a principal quantum number of 3 or 
greater. All these electrons and ions quickly thcrmalizc 
with the liquid helium. The ions form helium "snowballs" 
in a few picoseconds [g^ , and they do not move apprecia- 
bly from the sites where they are originated. On the other 
hand, as the energy of the free electrons drops below 
about 20 eV, the only process by which they can lose en- 
ergy is elastic scattering from helium atoms. Due to the 
low energy-transfer efficiency (about A = 1.36 x 10~^ per 
collision), these electrons make many collisions and un- 
dergo a random walk till their energy drops below 0.1 eV, 
the energy thought to be necessary for bubble state for- 
mation. Once thermalized, the electrons form bubbles in 
the liquid typically within 4 ps [zO] ■ Due to the Coulomb 
attraction, electron bubbles and helium ion snowballs re- 
combine in a very short time and lead to the production 
of HeS excimer molecules 



the charge extraction from a beta particle track, simu- 
lated by Guo et al. (t^ . and that from an alpha particle 
track, simulated by Ito et al. [ill, are shown as the blue 
solid curve and the red dashed curve, respectively. Note 
that in the low field regime, the measured charge col- 
lection by Ghosh [s^ and Sethumadhavan 81 1 for beta 
particles is higher than the predicted result by Guo et 
al. [T^I- Furthermore, these charge extraction analyses 
are for temperatures above 1 K. At very low tempera- 
tures, the ionized electrons can stray away from the field 
lines which enhances the charge extraction at a given ap- 
plied field [zll. 
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When an external electric field is applied, some of the 
electrons can escape the recombination and be extracted. 

At temperatures above 1 K, electron bubbles essen- 
tially move along the electric field lines in the moving 
frame of the ions due to the viscous damping [7l|, [t^ . 
In this situation, the fraction q of the electrons that can 
be extracted under an applied field £ depends largely on 
the initial electron-ion separation and the ionization den- 
sity along the projectile track. The mean electron-ion 
separation has been determined to be about 60 nm for 
both beta particle ionization events a-nd alpha par- 
ticle ionization events [73| . The energy deposition rate 
for an electron of several hundred keV is approximately 
50 eV/micron, whereas for an alpha particle of a few 
MeV the rate is 25 keV/micron [TJ]. The average energy 
needed to produce an electron-ion pair has been mea- 
sured to be about 42.3 eV for a beta particle [zll and 
about 43.3 eV [t^ for an alpha particle. It follows that 
charge pairs are separated on average about 850 nm along 
a beta particle track and only about 1.7 nm along the 
track of an energetic alpha particle. The recombination 
along a beta particle track where the electron-ion pairs 
are spatially separated is described by Onsager's gemi- 
nate recombination theory [t^I- For the highly ionizing 
track of an alpha particle in liquid helium, the electrons 
feel the attraction from all nearby ions and are harder to 
be extracted. Jaffe's columnar theory of recombination 
is more applicable in this situation [7^, [t^. In Fig. [51 



FIG. 8: (color online). Electron extraction fraction g as a 
function of applied electric field. The blue solid curve repre- 
sents the simulated electron extraction from beta tracks by 
Guo et al. [t^]. The red dashed curved represents the simu- 
lated electron extraction from alpha tracks by Ito et al. [t^. 

The mean electron-ion separation along the track of a 
low energy recoil He nuclei should be similar to that for 
beta and alpha particles. The ionization density along 
the He nuclei track can be estimated by {N^i + ^Nex)/Z, 

where Z= dE' /S{E') is the track length of the recoil 
He ion. Due to the Lindhard effect, a major part of the 
projectile energy is lost to elastic collisions at small E. 
Consequently the ionization density along the track of a 
recoil He ion should be much lower than that along the 
track of an energetic alpha particle. For instance, for a 
10 keV recoil He nuclei, the ions produced are on aver- 
age separated by about 20 nm along the track. At lower 
recoil energies, the separation between ionizations be- 
comes comparable or even larger than the mean electron- 
ion separation. As a consequence, the charge extraction 
fraction q for low energy nuclear recoils is expected to 
be similar to that for electron recoils. Due to the lack 
of experimental information, in the following analysis we 
assume the same q for both low energy recoil He nucleus 
and beta particles. The charge counts S2 for nuclear re- 
coils and electron recoils are thus given by q{Yei + ^Yex)E 

and q^Y^f' -\- ^Yex')E.^ respectively. Note that the terms 

^ ^e-L account for the ionizations produced by 
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the auto-ionization of the excited He atoms. 



2. Excitations and scintillation 

Excited hehum atoms can be produced in excitation 
coUisions. For electron recoils, Sato et al. [i^ calculated 
that 83% of the excited helium atoms produced in ex- 
citation collisions are in the spin-singlet states and the 
rest 17% are in triplet states. For low energy nuclear 
recoils, however, the direct excitations are nearly all to 
spin-singlet states [H, HI] . This is because that since the 
total spin is conserved, excitation to triplet states can oc- 
cur only when both the recoil He and the ground state He 
atom are excited simultaneously to triplet states. This 
process requires more energy and has a lower chance to 
occur. The excited atoms are then quickly quenched to 
their lowest energy singlet and triplet states, He*(2^S) 
and He*(2'^S), and react with the ground state helium 
atoms of the liquid, forming excited Y{c2{A^Y,u) and 
He2(a'^I]ii) molecules 



He* -t- He 



Hen 



(14) 



Hej excimer molecules are also produced as a con- 
sequence of recombinations of electron-ion pairs. For 
geminate recombination, experiments [zl] indicate that 
roughly 50% of the excimers that form on recombina- 
tion are in excited spin-singlet states and 50% are in 
spin-triplet states. He2 molecules in highly excited sin- 
glet states can rapidly cascade to the first-excited state, 
He2(A^I]„), and from there radiatively decay in less than 
10 ns to the ground state [s^, He2(X^Sg), emitting ul- 
traviolet photons in a band from 13 to 20 eV and centered 
at 16 eV. As a consequence, an intense prompt pulse of 
extreme-ultraviolet (EUV) scintillation light is released 
following an ionizing radiation event. These photons can 
pass through bulk helium and be detected since there is 
no absorption in helium below 20.6 eV. 

The radiative decay of the triplet molecules He2(a'^I]ti) 
to the singlet ground state He2(^^Sg) is forbidden since 
the transition involves a spin flip. The radiative lifetime 
of an isolated triplet molecule He2(a'^I]„)has been mea- 
sured in liquid helium to be around 13 s [83[. The triplet 



molecules, resulting from both electron-ion recombina- 
tion and from reaction of excited triplet atoms, diffuse 
out of the ionization track. They may radiatively decay, 
react with each other via bimolecular Penning ioniza- 
tion [U, or be quenched at the container walls. Experi- 
mentally, these molecules can be driven by a heat current 
to quench on a metal detector surface and be detected as 
a charge signal [H, [s^ . 

Note that non-radiative destruction of singlet excimers 
by the bimolecular Penning ionization process can lead 
to the quenching of the prompt scintillation light. Such 
quenching effect has been observed for energetic alpha 
particles in liquid helium [t^, H^j • At temperatures above 
2 K, the singlet excimers can diffuse on the order of 10 
nm in their 10 ns lifetime [881. Based on the discussion 



presented in the previous section, the mean separation of 
the excimers along the track of a low energy recoil helium 
can be greater than the diffusion range of the singlet 
excimers. The quenching of the prompt scintillation for 
low energy nuclear recoils should thus be small. At low 
temperatures, the quenching effect may be significant. 
However, it has been observed that even for the highly 
ionizing track of an energetic alpha particle, the light 
quenching becomes mild below about 0.6 K [13] • This 
is presumably due to the trapping of the excimers on 
quantized vortex lines that are created accompanying the 
energy deposition of the recoil helium [s^ . Such trapping 
limits the motion of the excimers and hence reduces the 
light quenching. Lack of experimental knowledge about 
the decay rates at which bimolecular Penning processes 
occur among the different excimers, we shall not include 
the quenching effect in our analysis. 

Based on the above discussion, the prompt scintillation 
photons (51 ('^^) and triplet molecules (S'3('^^) produced 
by an electron recoil event are given by 



S-l^") =E- 



Yef^ ■ (1 - 9) • 50% + • 86% • - 



+Y^^ ■ 86% ■ i • (1 - 9) 



(15) 



Yel^ • (1 - g) • 50% + ■ 14% • 



+Y^^) .U%■{l-q)■]- 



m 



The factor 2/3 accounts for the fraction of the excited 
atoms that do not undergo autoionization. The above 
two formulas assume that the recombination of electron- 
ion pairs produced by the autoionization of singlet (or 
triplet) helium atoms tends to generate only singlet (or 
triplet) helium excimers. The justification for this as- 
sumption is that the energy of the electrons produced 
in the Hornbeck-Molnar process is low (less than 2 eV). 
These electrons do not move very far from their parent 
ions, hence their spin correlation with their parent ions 
is likely strong enough to survive the whole ionization- 
recombination process. As for the nuclear recoils, the 
51*^"^ and 5'3("^ counts are given by 



5*1^") =E ■ 



Y,i • (1 - g) • 50% + y,,. • 3 



• 3 • (1 - <z) 



53^") = E-Y,i • (1 - <?) • 50% 



(17) 



(18) 



Since the excited atoms are assumed to be all in singlet 
states for nuclear recoils, the triplet molecules are gen- 
erated solely as a consequence of the recombination of 
charge pairs produced in direct ionization collisions. 
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For the readers' convenience, in tabic U we summarize 
the formulas that we used to estimate the SI, S2, and S3 
counts for both nuclear recoils and electron recoils with 
incident energy of E. 

C. Discrimination of nuclear recoil and electronic 
recoil 

1. Ratios of the signals from different channels 

The success of a direct dark matter experiment will 
depend in its ability to distinguish between electron re- 
coils and nuclear recoils. Discrimination between both 
types of recoils can be done by looking at the ratio of 
the counts from different signal channels. These ratios 
depend on the event type, the recoil energy, and the ap- 
plied electric field. The formulas listed in table U allow 
us to estimate these ratios. As an example, in Fig. |9l 
the calculated ratios of S2/S1 and S3/S1 are shown as a 
function of the applied electric field for both the electron 
recoils and nuclear recoils with a recoil energy of 10 kcV. 
The S2/S1 ratio for both electron recoils and nuclear re- 
coils increases with the applied electric field. This is be- 
cause at higher fields more electrons can be extracted, 
which enhances the S2 counts and at the meanwhile re- 
duces the SI counts since less electrons recombine with 
the ions to form singlet molecules. The difference of the 
S2/S1 ratio between electron recoils and nuclear recoils 
becomes greater at higher fields, which means that better 
discrimination based on S2/S1 can be achieved at higher 
fields. 
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FIG. 9: (color online) The ratio of the counts from different 
signal channels for 10 keV nuclear recoil and electronic recoil 
events as a function of the applied electric field. 

In Fig. [TOl we show the calculated ratios of S2/S1 and 
S3/S1 as a function of the event energy for both electron 
recoils and nuclear recoils at an applied field of 8 kV/cm. 
Since we take the ionization and excitation yields for elec- 
tron recoils to be constants, the S2/S1 and S3/S1 ratios 
for electrons recoils are independent of the recoil energy. 



0.8-1 ■ 1 ■ 1 • 1 • 1 • 1 

Drift Field: 8000 V/cm 

0.6- 

■| 0.4- ~ 

CO 

q: 

S2/S1 S3/S1 

0.2- ' .•• Nuclear recoil: — 

Electron recoil: 

0.0-1 ■ 1 ■ 1 ■ 1 • , • 

20 40 60 80 100 

Energy (keV) 

FIG. 10: (color online) The ratio of the counts between dif- 
ferent signal channels for nuclear recoil and electronic recoil 
events as a function of the event energy. The applied electric 
field is 8 kV/cm. 

For nuclear recoils, both the S2/S1 and S3/S2 ratios de- 
crease with decreasing recoil energy. Note that the ratios 
evaluated here are based on the calculated average counts 
from the different signal channels. In real experiment, 
there always exist number uncertainties of the counts. 
At low recoil energies where the counts are small, the 
relative uncertainties of the counts as well as the ratios 
of the counts between different channels become large, 
which limits the discrimination of the two types of re- 
coils. For helium detector, as we can see from Fig. [101 
the S2/S1 and S3/S1 curves for nuclear recoils bend away 
from those for electron recoils, which compensates the 
effect due to count uncertainty. As we shall show later, 
excellent event discrimination can still be achieved even 
down to a few keV energy regime. 

2. Scintillation efficiency factor 

The quantities that can be measured experimentally 
for a recoil event are the counts from the different sig- 
nal channels. One can plot, for instance, the S2/S1 ratio 
against the SI counts. However, the conversion between 
SI counts to the event energy is different for electron re- 
coils and nuclear recoils. For electron recoils, the event 
energy is proportional to the SI counts, since the ioniza- 
tion and the excitation yields are taken to be constant. 
For nuclear recoils, such conversion is nonlinear. The 
effective scintillation efficiency L^ff describes the differ- 
ence between the amount of energy measured in the de- 
tector between both types of recoils. In the notation used 
in the field, the keV electron equivalent scale (keVee) is 
used to quantify a measured signal in terms of the energy 
of an electron recoil that would be required to generate it. 
Similarly the keVr scale is used for nuclear recoil events. 
For a nuclear recoil of energy Er , the electron recoil event 
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TABLE I: The yields of prompt scintillation (SI), charge (S2), and He2 triplet molecules (S3) for nuclear recoils and electron 
recoils with incident energy of E in liquid helium. 





Nuclear recoils 


Electron recoils 


SI 


E ■ [O.S-Yei • (1 - g) + 0.67- + O.SS-Y;. ■ (1 - g)] 




• (1 - g) + + 0.29-y4"^ 


(1- 


g)l 


S2 


£-(n,+0.33-ye.)-q 


i5-(yi,^'+0.33-ye^')-g 


S3 


■ (1 - g) 




(1 - g) + 0.093- Yi,"' + 0.047- Yil' 


■(1 


-g)] 



that would produce an equivalent SI signal is given by 



Ee{keVee) = Lef fx.Er{keVr). 



(19) 



By definition L^f / is the nuclear recoil scintillation ef- 
ficiency relative to that of an electron recoil of the same 
energy at zero field. Experimentally, the conversion be- 
tween SI and the electron equivalent scale keVee can be 
established using gamma line sources, for example the 
^''Co 122 keV line. The nuclear recoil response as a func- 
tion of energy can be established using neutron scatter- 
ing, either in a mono-energetic neutron scattering exper- 
iment, or by using a neutron source with a broad energy 
distribution and comparing the observed shape of the 
nuclear recoil spectrum with detailed Monte Carlo sim- 
ulations. Using the formulas listed in Table U we can 
estimate the L^f / by calculating the ratio of the energies 
of the two types of recoil events that give the same SI 
counts. The result is shown in Fig. [TTJ 
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FIG. 11: (color online) The effective quenching factor Leff 
as a function of the recoil event energy. The x represents the 
calculated ie// for helium under zero applied electric field. 
The measured data for liquid Xenon by G. Plante et al. [s^ 
(A) and by A. Manzur et al. [9OI {+) are also shown. 

The event discrimination ability of a detector drops 
drastically below a certain threshold SI counts. For a 
given energy threshold in keVee scale, a detector with 
a higher L^ff has lower nuclear recoil energy thresh- 
old, hence would be sensitive to low energy WIMPs. In 
Fig. 1111 we also show the experimentally measured ^e// 
data for liquid Xenon [sol |90 | . In the low energy regime 
of a few keV, Xenon-based detector only has a L^ff of 



less than 0.1 while helium detector has a L^ff above 0.4. 
So while liquid helium has substantially lower scintilla- 
tion yield for electron recoils, this is unlikely to be the 
case for nuclear recoils. 



3. Rejection power 

The uncertainty of the signal counts limits the dis- 
crimination between the nuclear recoils and the electron 
recoils at low energies. To study this effect, we performed 
a Monte Carlo simulation. For a given recoil energy E 
in electron equivalent energy scale, we randomly gener- 
ate SI and S2 counts for a nuclear recoil event and an 
electron recoil event, assuming a Poisson distribution of 
the counts with mean count values given by the formulas 
listed in Table HI In each trial, the ratios of S2/S1 for a 
nuclear recoil and for an electronic recoil are evaluated 
and represented by a red dot and a black dot in the S2/S1 
versus E plot. 10^ trials are carried out for each energy. 
An example is shown in Fig. 1121 To match with real ex- 
periments, we assume that only 20% of the scintillation 
photons (SI) are collected (typical for a two- phase detec- 
tor as we shall discuss later), and that all the extracted 
electrons (S2) under the applied drift electric field can 
be detected. A clear separation can be seen between the 
electron recoil band and the nuclear recoil band, a neces- 
sary criterion for any direct dark matter experiment. At 
low energies, the two bands overlap due to the large scat- 
tering of the S2/S1 value. This large scattering is caused 
by the large relative uncertainty of the counts when their 
mean values of the Poisson distributions are small. 

To calculate the rejection power, we divide the two 
bands in energy slices. We select the lower half of the 
nuclear recoil band as our WIMP region of interest. At 
low energies it is crucial to know what percentage of the 
electronic recoil band leak into the lower 50% nuclear re- 
coil band. The rejection power (or sometimes called dis- 
crimination power) is found as the fraction of electronic 
recoil events below the nuclear recoil centroid. A full 
description on this method has been given by A. Man- 
alaysay |9l|. In Fig [131 we show the calculated rejection 
power as a function of event energy in keVr scale at sev- 
eral applied electric fields and with different SI collection 
efficiencies. At low energies, the ability to distinguish 
electron and nuclear recoils is degraded because of the 
lack of charge and light signal. But above a few keV, dis- 
crimination power is predicted to improve considerably, 
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FIG. 12: (color online) Monte Carlo simulation of S1/S2 ra- 
tio for nuclear recoils (red dots) and electronic recoils (black 
dots). The SI scintillation light collection efficiency is as- 
sumed to be 20%. The applied electric field is 8 kV/cm. The 
event energy is in the keV electron equivalent energy scale 
(keVee). 



and this should allow for low-background operation and 
a sensitive WIMP search. The discrimination is better at 
higher fields or with higher SI collection efficiency. We 
considered fields up to 40 kV/cm. It has been shown 
that such high fields can be readily applied in liquid he- 
lium [t^. Indeed, the design electric field value of the 
Spallation Neutron Source (SNS) neutron EDM experi- 
ment is 50 kV/cm (92| - [93 |. As we shall discuss later, for 
a single-phase helium detector operated at very low tem- 
peratures, sensitive bolometers immersed in liquid helium 
may be used to read out the light and the charge signals. 
In this case, 80% SI collection may be possible with the 
detector inner surface being covered by bolometer arrays. 
Note that the rejection power analysis is based on the 
the charge extraction curve shown in Fig. [51 The actual 
charge extraction at a given field could be higher, espe- 
cially at low temperatures where the ionized electrons 
can stray away from field lines [tI] . Considering this fac- 
tor, the actual rejection power could be better than the 
results shown in Fig |13l 



III. HELIUM-BASED DETECTOR 

For dark matter detection, we propose to detect both 
prompt scintillation and charge in liquid helium-4, us- 
ing a time projection chamber design. This is essent ially 
the same approach used in Ar and Xe detectors j95l - l97j , 
which has proven to be very effective, providing excellent 
position resolution and gamma ray discrimination. Based 
on the readout schemes for the light and charge signals, 
we discuss two proposals for a liquid-hclium-based dark 
matter detector. One proposal is to operate the dctcc- 
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FIG. 13: (color online) Calculated rejection power for a he- 
lium detector as a function of event energy in keVr scale. 



tor at high temperature regime 3 K) using photo- 
multiplier tube arrays for signal readout, and the second 
proposal is to operate the detector at low temperatures 
(~ 100 mK) using bolometer arrays for signal readout. 



A. High temperature operation scheme 

Operating a helium detector at relatively high temper- 
atures is favored in economy since such a detector does 
not require complicated large-scale dilution refrigeration 
system. At high temperatures where the helium vapor 
density is high, some existing technologies for charge sig- 
nal amplification may be ready applied to the helium 
detector, such as the proportional scintillation in a two- 
phase chamber that has been used in Argon, Krypton, 
and Xenon detectors [95l - [98| . or the ele ctron avalanche 
in a Gas Electron Muhiplier (GEM) f99l-[l01|. 



A conceptual schematic of a two-phase helium-based 
time projection chamber is shown in Fig. 1141 Ionizing 
radiation events in liquid helium produce both prompt 
scintillation light (SI) and ionizations. A drift electric 
field can be maintained between the anode and the cath- 
ode. Some ionized electrons can be extracted into the gas 
phase and caused to emit 16 eV EUV photons (S2) by a 
strong field in the vapor. The anode and cathode may 
be a transparent material coated with Indium Tin Oxide 
such a s rece ntly demonstrated in the DarkSide-10 exper- 
iment [l02| . so produce a uniform proportional scintil- 
lation field and eliminate liquid helium scattering events 
below the cathode. To detect the EUV photons, as is typ- 
ical for scintillation detection in liquid argon, liquid neon, 
and liquid helium, the inner surface of the chamber win- 
dows and the transparent electrodes can be coated with 
tctraphcnyl butadiene (TPB) wavelength shifter. The 
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FIG. 14: (color online). A schematic of a two-phase helium- 
based dark matter detector. 



EUV scintillation light is converted to the visible on the 
TPB coating, which has approximately 100% photon-to- 
photon conversion efficiency at the helium scintillation 
wavelength of 80 nm [H. Photomultiplier tube (PMT) 
array placed outside the helium chamber can be used 
to detect the converted photons. SI light detection ef- 
ficiency of about 20% may be expected in such a detec- 
tor, sim ilar to that measured by the DarkSide collab- 
oration 



102[, which recently demonstrated a zero- field 
signal yield of 9.1 photoelectrons/keVee in a two-phase 
liquid argon detector. With time projection readout, the 
time between the SI and S2 signals indicates the depth 
(z) of the event, while the hit pattern in the top array 
of photomultipliers gives the x — y position of the event. 
This allows determination of a low-background central 
fiducial volume, which is used for the dark matter search. 
Events close to chamber surface may be discarded, and 
the S2/S1 ratio provides discrimination power against 
gamma ray background. Any gamma ray or neutron 
events that cause multiple scatters will generate multi- 
ple S2 signals, and these may also be discarded. 

The extracted electrons moving in helium vapor un- 
dergo collisions with helium atoms. Due to their very 
small mass the electrons can give up almost no energy 
to the helium atoms in the course of elastic collisions. 
When the kinetic energy of the electrons builds up over 
a few mean free pathes to exceed the excitation thresh- 
old 20 eV) of helium atoms, inelastic collisions be- 
tween the electrons and the helium atoms, which lead 
to the production of excited helium atoms, can occur. 
The excited helium atoms in spin-singlet states can react 
with surrounding helium atoms and decay to the ground 
state by emitting 16 eV EUV photons (S2). The S2 
strength increases with the applied field. However, when 
the electric field in vapor is too strong, the electron en- 
ergy can exceed the ionization threshold 24.6 eV) of 



helium atoms. In this situation, charge multiplication 
in gas occurs, and eventually avalanche breakdown can 
take place. Given the premium on high electric field for 
getting good event discrimination, it is advantageous to 
operate the detector at a field in the vapor only slightly 
below the breakdown field. The breakdown voltage Vbd of 
helium gas in a uniform field generated by a pair of elec- 
trod es se parated by a distance d is given by the Paschen's 
law fl03 



Vm = 



A-p-d 



ln{p-d) + B 



(20) 



where p is gas density. A and B are experimentally de- 
termined constants. Some available experimental data 
of the brea kdown vo ltage for helium gas are shown in 
Fig. [15] (a) |l04 Il05{ . The solid curve represents a typi- 
cal Paschen's curve for helium gas obtained by fitting the 
experimental data using Eg. fIU\ 105, 106]. 
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FIG. 15: (color online), (a) The dielectric breakdown volt- 
age in helium g function of the product of the gas 
density p and the electrode separation d. The red circles 
(o) and t he black squa res (□) are experimental data from 
Ref. [l04l ] and Ref. [lOSi ]. respectively. The solid curve repre- 
sents the Pasch en's curve obtained by fitting the experimental 
data [lOSl . Il0(j | . (b) The dielectric breakdown field for satu- 
rated helium vapor as a function of temperature. 
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For saturated helium vapor in a two-phase hehum de- 
tector, the gas density as a function of temperature is 
known |l07l |. We can thus use the Paschen's curve to 
calculate the breakdown field for a two-phase helium de- 
tector as a function of the operation temperature for a 
given electrode separation d in the vapor. In Fig. [15] (b), 
the calculated breakdown field at d = 5 mm is shown. 
As wc can sec, the breakdown field increases drastically 
with temperature due to the increased vapor density. At 
3.2 K, the breakdown field is about 40 kV/cm. Under 
a drift field i^drift that is close to this breakdown field, 
good event discrimination power is expected. 

The drift speed of the electrons v in liquid helium is 
given hy V = fJ-Edrut, where fi is the electron mobility in 
the liquid. Electrons drift much slower in liquid helium 
than in liquid xenon and liquid argon due to the small 
/i associated with the bubble state in helium. Ba sed on 
the known mobility of the electrons in helium |l08j , their 
drift velocity under a field close to the vapor breakdown 
field as given in Fig. [13] (b) can be calculated. The result 
is shown in Fig. 1161 The electron speed drops with in- 
creasing temperature below the lambda point (~ 2.17 K), 
which is due to the steep drop of the electron mobility. 
Above the lambda point, fj, decreases slowly with increas- 
ing temperature, and the electron speed rises steadily as 
the breakdown field increases. At 3.2 K, the electron 
speed is about 10 m/s. For a 10-liter sensitive volume, 
if we take the distance from the bottom electrode to the 
liquid/gas interface to be 20 cm, the maximum delay 
time between SI and S2 will be about 20 ms. For a well 
shielded detector, event pileup should be well below the 
level that would cause mismatching of SI and S2 signals. 



electrons using GEMs has already been s tudied ex per- 
imentally by the "e-bubble" cohaboration [lolini- A 
disadvantage of the GEM is that it gives poorer energy 
resolution than proportional scintillation since it relies 
on a breakdown for electron gain. However, because of 
the slow electron drift speed, the electrons will arrive at 
the GEM one at a time, easily distinguishable due to the 
excellent GEM timing resolution. We may operate the 
GEM in a single electron detection mode, counting single 
electron pulses instead of using the pulse sizes to deter- 
mine the event energy. Note that it was shown that in ul- 
trapure heliu m gas GEM can operate only at charge gain 
close to unity jllll . Ill2l |. However, during the avalanche 
development, excited helium atoms and molecules are 
produced. The decay of these electronic excitations leads 
to the emission of 16 eV photons. Instead of detecting 
the charge produced in the GEM, we may detect these 
photons with arrays of PMTs. At the same time the 
GEM (or stack of GEMs) could also be used to amplify 
the prompt scintillation signal. The side of the GEM 
facing the liquid could be coated with Gesium Iodide 
(Gsl) |ll3l | or other photocathode material so as to be 
sensitive to the prompt scintillation light. Furthermore, 
the extracted electrons may also be detected via electro- 
luminescence produced under very high (^1-10 MV/cm) 
fields on the surface of thin wires or a GEM immersed 
in liquid helium. Such ele ctrol uminescence has already 
been observed in liquid Ar 114 1 and liquid Xe [l3l ■ This 
could allow electrons to be individually detected, while 
not subjecting gaseous helium to a strong electric field. 



B. Low temperature operation scheme 



17.5^ 




2.0 2.5 3.0 3.5 4.0 
T(K) 

FIG. 16: (color online). The electron drift velocity in liquid 
helium as a function of the temperature under a drift field that 
is close to the breakdown field for saturated helium vapor. 

Alternatively one may detect the extracted electrons 
using Gas Electron Multipliers (GE Ms) o r Thick Gas 
Electron Muhipliers (THGEMs) fool-floH. Detecting 



At very low temperatures (e.g. 100 mK), calorimetric 
sensors with small heat capacity can be used for par- 
ticle/photon detection with remarkable sensitivity and 
low threshold. A significant advantage of using calori- 
metric sensors is that one could in principle cover all the 
inner surface of the detector with calorimetric sensor ar- 
rays, while not being limited by the 20-30% quantum 
efficiency typical of photocathodes. SI light collection 
efficiency approaching 100% might be achieved, which 
would improve the rejection power of the detector. At 
low helium temperatures, the thermal coupling between 
the calorimetric sensors and the liquid helium is weak, 
enabling them to be immersed in the liquid helium with- 
out losing much thermal signal to the bath. 

The field of low temperature bolometers is develop- 
ing rapidly. The types of temperature sensors most 
commonly used are neutron-transmutation doped (NTD) 
Germanium thermistors, and superconducting transition 
edge sensors (TES) till [nil- The electric conductivity 
of NTD sensors strongly depends upon the temperature, 
with typically resistance of 1 '--^ 100 Mfl at low tempera- 
tures. NTD thermistors are easy to use because they can 
be operated with conventional electronics. Mass produc- 
tion is also possible for the NTD sensors. A TES is a 
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superconducting strip operating at the temperature of 
its superconducting-normal transition. The resistance in 
the normal state is usually 10 mfi ~ lil, and the temper- 
ature dependence of the resistance can be very large at 
the transition. Recent developments include not only im- 
provements in single TES's, but also large arrays of TES's 
and techniques for multiplexing them. For this dark mat- 
ter application, especially promising are non-equilibrium 
detectors, in which interactions produce quasiparticles 
in a superconducting st rip, which can be then collected 
in a TES and detected [ll7l | . Another possible choice is 
metallic magnetic calorimeters (MMC) [lT8l4l20t . MMCs 
are based upon the use of magnetic sensors to measure 
very small temperature changes resulting from the ab- 
sorption of energy by energetic photons. Instead of mea- 
suring the resistance of a sensor such as an NTD or a 
TES, an MMC measures the change of magnet ization of 
paramagnetic ions in a metallic host |l2ll . Il22| . The in- 
ternal thermalization time of MMCs is very fast, within 
a microsecond. 
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FIG. 17: (color online). A schematic of a single- phase helium 
detector operated at mK temperatures. 



A conceptual schematic of a calorimeter-based single- 
phase helium time projection chamber is shown in 
Fig. [171 The prompt 16 eV photons produced by a recoil 
event hit the bolometer arrays and deposit their energy 
in the sensors which give rise to the SI signal. At low 
temperatures, bolometers may be made sensit ive enough 
to allow the detection of individual photons [llTl ll2C| . 
To extract the ionized electrons, three electrodes are ar- 
ranged in a way to drive the electrons toward either the 
top or the bottom electrode. Since all the electrodes arc 
immersed in liquid helium, high voltages ±100 kV) 
can be applied to them. The drift field with the shown 
electrode arrangement can be made as high as a few tens 
of kV/cm, which helps to improve the rejection power 
of the detector. To amplify and detect the charge sig- 
nal, both the top and the bottom electrodes are made 
of thin wire arrays. The extracted electrons drifting to- 
wards these thin wires may produce electroluminescence 



as they approach the wire surface. These photons (S2) 
can be detected by the same bolometer arrays. Again, 
event position reconstruction can be made based on the 
delay between SI and S2, and the hit pattern on the top 
(or bottom) bolometer array. 

Note that the mobility of electron bubbles in liquid he- 
lium increases drastically with decreasing temperature. 
Under saturated vapor pressure, if the velocity of the 
electrons exceeds a threshold velocity of the order of 40 
m/s, quantized vortex rings are nucleated. The electrons 
can get stuck on them, and the charged vortex ring moves 
through the liquid as a single entity. A striking feature 
of the electron-ring con iplex is that its velocity decreases 
with increasing energy |123l | . When a strong drift field (a 
few kV/cm) is applied, the velocity of the cha rged vortex 



rings can be as low as on the order of 1 m/s |124{ . How 



ever, it has been shown that at low temperatures when 
isotropically purified helium is pressurized to above 15 
bar, electrons can be driven at a speed close to or higher 
than the Landau velocity (~ 50 m/s). Instead of nucleat- 
ing vortex rings, the electrons spontaneously emit roton 
pairs [l25l - [r27ij . The rate of roton emission depends on 
the field strength. Furthermore, it has been shown that 
when the electron speed is not too much higher than the 
Landau velocity, the majority of the emitted rotons tend 
to have momentu m ah gned in the same direction with the 
electron velocity |l28j . A roton beam is formed accompa- 
nying every extracted electrons. Note that rotons in the 
R"'' branch move along the electron velocity direction due 
to their positive effective mass, while R~ rotons are emit- 
ted in the opposite direction since their effective mass is 
negative in helium jl29j . Despite the low transmission 
coefficient of the roton energy across the liquid-solid in- 
terface at low temperatures, a fraction of the roton en- 
ergy can ne verth eless transmit into the bolometers and 
be detected [l30j | . Detecting the rotons provides a unique 
way for charge signal amplification and detection, with 
potentially hundreds of keV of roton signal produced by 
each drifted electron. Operation at pressure > 15 bar, as 
required in this approach, may also be advantageous for 
maintaining higher drift field by suppressing gas bubble 
formation. 

Note that in the initial proposal by Lanou et al. [soj . 
the idea of charge signal amplification via roton/ vortex 
generation was briefly mentioned. It was proposed that 
event discrimination might be achieved based on the ratio 
of prompt rotons accompanying the initial recoil deposi- 
tion to the delayed rotons from the drifted charge. Due to 
the low transmission of roton energy into the bolometer 
surface, detecting the prompt rotons for low energy re- 
coils can be quite challenging. Nevertheless, if detection 
of prompt rotons and phonons could be accomplished 
for very low energy nuclear recoil events, search for ex- 
tremely light WIMPs may be conducted. At very low 
energies, ionization is strongly suppressed for nuclear re- 
coil events, and even a single electron from an electron 
recoil event would produce a large roton signal, allowing 
electron recoil backgrounds to be vetoed. In addition. 
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the roton/phonon ratio may be different for nuclear and 
electron recoils, allowing electron recoil backgrounds to 
be discriminated by bolometer pulse shape. 

It is worthwhile mentioning that at low temperatures, 
metastable helium molecules in triplet state can drift at 
a speed of a few meters per second [11]. When these 
molecules collide on the bolometer surface, they undergo 
non-radiative quenching and release over 10 eV of energy 
depending on the material of the impinged surface. A 
significant amount of this energy will be deposited into 
the calorimctric sensor, which may allow us to detect the 
molecule signal (S3). A combined analysis of S2/S1 and 
S3/S1 ratios may further improve the rejection power of 
the detector. 



IV. SENSITIVITY 

In any direct dark matter detection scheme, the pri- 
mary requirement is the strong reduction of radioactive 
backgrounds. The approaches described above are de- 
signed to enable this, since liquid helium may readily be 
purified of impurities, and because the ratios of signal 
channels may be used to identify gammas and betas on 
an cvent-by-event basis. 

Internal backgrounds (due to radioactive impurities 
within the target material) are particularly straightfor- 
ward to eliminate in liquid helium. Like other noble 
gases, helium is readily purified with heated getters to 
remove any chemically reactive species, which includes 
anything that is not a noble gas. In addition, helium has 
no long-lived radioactive isotopes and therefore has no in- 
trinsic backgrounds, unlike argon and krypton which con- 
tain the beta emitters ^^Ar and ^^Kr. Activated charcoal 
adsorber may then be used to remove all other radioac- 
tive noble gases (e.g. radioactive argon, krypton, and 
radon), since all other noble gases have larger polariz- 
abilities and masses, leading to larger binding energies to 
charcoal and substantially larger adsorption coefficients 
[l3lL Il32l | . Getters and charcoal may be used to purify 
the helium after it is transported underground. Any re- 
maining impurities will fall out of liquid helium and stick 
to plumbing and detector walls, and any beta, alpha, or 
nuclear recoil background due to these impurities may be 
removed through position reconstruction. 

External backgrounds may be reduced through shield- 
ing, careful detector materials selection, and event dis- 
crimination. The dominant background in this experi- 
ment is expected to arise from gamma rays that Comp- 
ton scatter in the helium and produce low-energy events 
that could be confused with dark matter particles. This 
gamma ray Compton scattering background tends to be 
flat with energy at low energies, and a typical Comp- 
ton scattering background rate for a shielded dark mat- 
ter experiment is about 1 event /keVee/ kg/ day, though 
this may be reduced further with special care given to 
shielding and detector materials. Neutron backgrounds 
are typically well subdominant to gamma rays, for most 



detector construction materials. Radon daughter back- 
grounds on inner detector surfaces can be troublesome, 
but are eliminated in this scheme through the excellent 
position reconstruction inherent to noble liquid time pro- 
jection chambers. 

With gamma rays generating the dominant back- 
ground, it is crucial to have excellent electron recoil ver- 
sus nuclear recoil discrimination. From the quantitative 
estimates described in Section II above, we have good 
confidence that liquid helium will indeed provide excel- 
lent discrimination power. In addition, the helium detec- 
tor may surrounded with a veto detector that is sensitive 
to gamma rays. Using organic scintillator, liquid xenon, 
or liquid argon as detector materials, such a veto may 
be used to efficiently tag gamma rays that small-angle- 
scatter in the helium fiducial volume, produce low energy 
events, and escape. Efficiencies of 90-98% may be ex- 
pected, as pr edicted fo r the DarkSide and LUX-ZEPLIN 
experiments [l02l Il33| . 




WIMP Mass [GeV/c^] 



FIG. 18: (color online). Spin-independent WIMP exclusion 
curves (solid lines), potential WIMP signals (solid regions), 
and projected liquid helium 90% sensitivity curves (dashed 
lines) in the region of 1-100 GeV WIMP mass. Exclusion 
curves include DAMIC in red [l35[, CDMS-II in green 
XENONIO in magenta [H, and XENONIOO in blue [H. Po- 
tential WIMP signals include DAMA in red [HI, CRESST 
in light blue \lT^ . and CoGeNT in green [l6| . Projected liq- 
uid helium sensitivities for different detector parameters are 
shown as dashed lines, including light blue (10 kV/cm, 20% 
SI collection, 4.8 keV threshold), green (20 kV/cm, 20% SI 
collection, 4.4 keV threshold), blue (40 kV/cm, 20% SI col- 
lection, 4.2 keV threshold), red (20 kV/cm, 80% SI collection, 
2.8 keV threshold), and magenta (40 kV/cm, 80% SI collec- 
tion, 2.6 keV threshold). Predicted limits assume an electron 
recoil background of 1 event/keVee/kg/day and a 95% effi- 
cient gamma ray veto. 
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As explained in section I, a useful figure of merit for 
light WIMP searches is (nuclear mass) x (energy thresh- 
old), which must be minimized to get the best light 
WIMP sensitivity. In the case of liquid helium, this must 
be balanced with the background reduction achieved 
through discrimination of electron recoil events, which 
improves with higher energy. Given helium's large pre- 
dicted nuclear recoil signals and excellent discrimination, 
we expect an energy threshold of about 4-5 keV with 
photomultiplier readout, potentially reducible to 1-3 keV 
with bolometric readout. The low nuclear mass of helium 
then gives access to very low WIMP masses, while still 
having significant background reduction through discrim- 
ination. 

While liquid helium will not provide significant self- 
shielding against gamma rays and neutrons (a sig- 
nificant background rejection method in LXe and 
LAr detectors), a plausible background rate of 10"'^ 
events/day /keVee/kg after discrimination will allow ex- 
cellent sensitivity to light WIMPs, for which current ex- 
perimental sensitivities are relatively weak. A detailed 
discussion of the background of a helium detector de- 
sig ned f or the HERON project was given by Huang et 
al. [l34l |. For a 1 kg helium fiducial mass, 20% light 
collection, a 20 kV/cm drift field, an energy threshold 
of 4.8 keV, 300 days of operation, and a 95% efficient 
gamma ray veto, one background event is predicted, with 
a WIMP-nucleon cross-section sensitivity of 10^"*^ cm^ at 
5 GeV, the dark matter mass predicted by asymmetric 
dark matter models. Sensitivity may be improved further 
with higher drift fields, more efficient light collection, and 



larger exposure, potentially reaching lO"^** cm^ or better 
between 2-20 GeV. Some predicted light WIMP sensitiv- 
ities are summarized above in Figure [THl 



V. CONCLUSION 

We conclude that liquid helium is an intriguing mate- 
rial for the direct detection of light WIMPs, as it com- 
bines multiple signal channels, comparatively large sig- 
nals for nuclear recoils, a low target mass, and the ca- 
pacity for electron recoil discrimination. In the detector 
schemes proposed here, a high electric field is used to 
extract electrons from nuclear recoil tracks, allowing a 
sizable charge signal, time projection chamber readout, 
and good position resolution. Before dark matter exper- 
iments can be performed with this technology, a method 
of detecting single electrons in liquid superfiuid helium 
must be demonstrated. In addition, detailed measure- 
ments must be done of the nuclear and electron recoil 
signal and discrimination efficiency at low energies. 
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